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ABSTRACT 

Observations on the high-redshift galaxies at z > 6 imply that their ionizing emissivity 
is unable to fully reionize the Universe at z ~ 6. Either a high escape fraction of 
ionizing photons from these galaxies or a large population of faint galaxies below the 
detection limit are required. However, these requirements are somewhat in tension 
with present observations. In this work, we explored the combined contribution of 
mini-quasars and stars to the reionization of cosmic hydrogen and helium. Our model 
is roughly consistent with: (1) the low escape fractions of ionizing photons from the 
observed galaxies, (2) the optical depth of Cosmic Microwave Background (CMB) 
measured by the WMAP-7, and (3) the redshift of the end of hydrogen and helium 
reionization at z ~ 6 and z ~ 3, respectively. Neither an extremely high escape 
fraction nor a large population of fainter galaxies is required in this scenario. In our 
most optimistic model, more than ~ 20% of the cosmic helium is reionized by z ~ 6, 
and the ionized fraction of cosmic helium rapidly climbs to more than 50% by z ^ 5. 
These results may imply that better measurements of helium reionization, especially at 
high redshifts, could be helpful in constraining the growth of intermediate-mass black 
holes (IMBHs) in the early Universe, which would shed some light on the puzzles 
concerning the formation of supermassive black holes (SMBHs). 

Key words: dark ages, reionization, first stars - intergalactic medium - quasars: 
general 


1 INTRODUCTION 

One major focus in the present cosmology is the reionization 
of the intergalactic medium (IGM) in the Universe. Obser¬ 
vations of the Cosmic Microwave Background (CMB) from 
the 7-yr Wilkinson Microwave Anisotropy Probe (WMAP- 
7) measure the Thomson electron scattering optical depth 
T = 0.088 ± 0.015, suggesting that hydrogen reionizatio n 
begins no later than z = 10.6 ± 1.2 jKomatsu et al.ll201ll) . 
Additional pieces of evidenc e come from the absorption spec¬ 
tra o f high-redshift quasar s (IFin et al.||2006l; iM ortlock et al.l 
I2OIII : iBolton et al.l I2OIII : ISchroeder L al.l l2013h and the 
rapid evo lution of the luminosity function of Ly a emitters 
at z > 6 dOuchi et al.ll201^ : ISchYker et al]l201^ . implying 
an end to hydrogen reionization at 6 < z < 7. All these con- 


* E-mail: haojm@mail.ustc.edu.cn; yfyuan@ustc.edu.cn; lei- 
wang@pmo.ac.cn 


straints suggest that hydrogen reionization is likely to be a 
much m ore extended process i n the redshift range 6 < z < 
12 (e.g. [Wyithe fc Loebll2003l: IChoudhurv fc Ferraral I2OO6I: 


lllev et al.l 20071 : Haardt &: Madaul 2012 : Robertson et al.l 

20131) 


Among various possible ionizing sources, stars in 
galaxies are commonly thought to be the most likely can¬ 
didat es (e.g. iFaucher-Ciguere et al.l I2OO8I : iRobertson et al.l 
I2OIOI) . However, the ionizing emissivity of the observed 
galaxies at z > 6 seems insufficient to maintain the 

cosmic hydrogen in a fully ionized state at z ~ 6 un¬ 
less the escape fraction of ionizing photons from these 
galaxies is larger than ~ 30%, otherwise a large popu¬ 
lation of fai nt galaxies below the detection limit would 
be required (iBolton fcjlaelmehj^ 2007 : Bunk er et al.| 2010|; 
Trenti e t al.l 2 010l: lAlv mez et^ 20 121: Bouwnnset_^lJ 


2^^~^ Fontanot et^20121: iFinkelstein et al 

Kuhlen fc Faucher-Ciguerd 20121 : Robertson et al.l 
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I Wise et al.ll20l4) . The requirement for such a high escape 
fraction at 2 > 6 is somewhat in tension with current obser- 
vational constraints on t he escape fra c tion a t lower redshifts 
llHeckman et alJ 2001: llwata et alJ 2 0091: Wnzella et alJ 


I 2 OI 0 I : iNestor et alJ 20111 " Vanzella et aLT2012[ l. Although 

there exists some evide nce supporting the redshift evolution 
of the escape fraction (iRazoumov fc Sommer-Larsen 200^ 


Paardekooper et al.ll201ll: Ferrara fc Loebll2013l : lMitra et alJ 


2 OI 3 I : ll^ntanot et alJ 2014ll . which is still inconclusive. On 

the other hand, some theoretical studies suggest that the 
star formati on activity in faint galaxies might be stron gly 
suppressed (iKrumholz fc Dekelll2()l iKuhlen et al.ir2012l l. 

Introducing the contribution from other ionizing sources 
might partially solve this so-called “photon-starved” prob¬ 
lem. Although it is widely believed that the contribution 
of quasars to the ionizing radiation decreases rapidly from 
2 ^ 3 to 2 ^ 6 and becomes negligible at 2 > 6 (e.g. 
iFan et ^l200ll : ICowie et al.ll2009l : IWillott et al.ll2010ll . the 
contribution from mini-quasars at high redshifts could still 
play an important role. The existence of black holes (BHs) as 
massive as ~ 10® M© at 2 ~ 7 becomes reasonable only when 
starting from a seed BH 0, which increases its mass via gas 
accretion or/and merging with each other. This scenario im¬ 
plies that there should be plenty of intermediate-mass black 
holes (IMBHs) with masses of the order of 10® —10® Mq shin¬ 
ing at 2 > 7. The recent identihcation of IMBHs also pro¬ 
vides strong evidence fo r this scenario llFarrell et al.ll20(^ : 
iLiitzgendorf et aP I 2 OI 3 I I . If such IMBHs accrete as mini¬ 
quasars at high redshifts, they could be important sources 
of th e early reionization of the IGM. A number of authors 
(e.g. iMadau et aDl2004 iDiikstra et ahlliool : IZaroubi et al.l 
|2007^ argued that mini-quasars, consistent with the con¬ 
straint from the soft X-ray background (SXRB), could still 
provide a relevant contribution to the reionization of the 
IGM at 2 > 6. 

Recently, by exploring the high-redshift Mbh — cr* rela¬ 
tion in light of the numerical simulation of the merger tree 
of dark m atter halos and us ing the Press-Schechter (PS) 
formalism, IWang et al.l (120101 1 calculated the central IMBH 
mass density as a function of redshift. And furthermore, by 
combining it with a parameterized UV photon emission ef- 
hciency, they found that the contribution of mini-quasars 
to hydrogen reionization lies in the range ~ 25% — 50% at 
2 ~ 6. However, they did not consider the impact of the 
helium ionizing radiation of these mini-quasars on the IGM 
reionization and ignored the contribution from stars. In this 
paper, we adopt the same BH mass density to estimate the 
number of ionizing photons from mini-quasars, assuming a 
power-law spectrum for the BH emission. In comb inatio n 
with the htting star formation rate density from Q ll2008ll . 
we assess the combined contribution of mini-quasars and 
stars to the reionization of cosmic hydrogen and helium. 

The structure of this paper is as follows. In Sect. 2, we 


^ The origin of the initial seed BHs remains an open question. 
The most obvious candida te is the collapse of massive Population 
III st ars at z > 20 (e.g. iM adau &: ReesI I 2 OOII : iHaiman fc LoebI 


ill st ars at z 2U (e.g. |Ma d 

l200ll : lTanaka fc Hairnm 2009l l~ 


, leaving behind remnant BHs with 
masses > 100 Mq. An alternative possibility is the direct col- 


> lO^K (e.g. 

Oh fc HaimanI l200d; iBromm fc LoebI l2003l; 

iLodato Sc Nataraian 

I 2 OO 6 I: iTanaka & Lill2014h. 


present the formalism for computing the reionization his¬ 
tories of hydrogen and helium and the number of ionizing 
photons from different types of sources. The implications of 
our results are discussed in Sect. 3, and our conclusions are 
summarized in Sect. 4. The cosmological p arameters used 
in th is paper are from the WMAP-7 results (iKomatsu et al.l 

I 2 OIIII : = 0 . 266 , Ha = 0.734, Hb = 0.0449, h = 0.71, 

erg = 0.801. We use n = 1 for a scale-invariant primordial 
power spectrum. 


2 THE REIONIZATION HISTORY OE 
HYDROGEN AND HELIUM 


The X-ray or UV emission from (mini-)quasars and stars can 
ionize the neutral hydrogen (H I), helium (He I) or singly ion¬ 
ized helium (Hell) in the IGM. In order to estimate the evo¬ 
lution of the regions of the ionized hydrogen and helium dur¬ 
ing the epoch of reionization, we follow the straightforward 
considerations of iBarkana fc Lo^ ll200lll . Given the number 
density of ionizing photons from stars and (mini-)quasars in 
a homogeneous but clumpy Universe, the evolution of the 
volume Hlling factor Qi (for species i = HH, Hell, Helll) of 
the ionized regions can be described by 


dQi 

dt 


1 

nf(l + 2)3 dt 


aB,iGne(l -I- e)^Qi, 


( 1 ) 


where n® is the comoving density, and = 

1.88 X 10”’'(Hbh®/0.022) cm"® and = 0.19 x 

10“^(Hbh®/0.022) cm"®, respectively. Here rie is the 
number density of electrons, as is the case B re¬ 
combination coefficient and C = (^h) /h-n the 
volume-averaged clumping factor of the IGM. The re¬ 
combination coefficients for HII, Hell and He III a t 
T 10^ K are adopted as follows dWvithe fc LoebIl2003l l: 
Q!B,Hii = 2.6 X 10"^® cm® s"®, QB.Heii = 2.73 X 10"®® cm® s"® 
and an,Belli = 13 x 10"®®cm®s"®, while for the clumping 
factor C, we use the following simple form given by 
iHaiman fc BrvarJ ll2006l') : C{z) = 1 -|- 9(7/(l -I- 2 ))® for 2 ^ 6 
and C — 10 for 2 < 6. The number densities of ionizing 
photons for HH, Hell and He HI are calculated through 
the integration of the corresponding spectra in the energy 
intervals 13.6 < E < 24.6 eV, 24.6 < E < 54.4 eV and 
E > 54.4 eV, respectively. Two classes of sources, mini¬ 
quasars and stars, are considered as the contributors to the 
ionizing photons, which will be discussed in the following 
sectio ns. We also note that, following IWvithe fc LoebI 
ll2003l) , two constraints are considered in this model, 
namely, (1) (Jneiii ^ Qhii; and (2) Qaeii ^ Qhii — QHeiii, 
showing that the helium ionization front never overtakes 
the hydrogen ionization front. 

The optical depth of CMB photons due to the Thomson 
Scattering is given by 


Te = y d2^C<TTne(2)(l -I-2)®, 


( 2 ) 


where (tt = 6.65 x 10 ®® cm ® is the cross section of Thom¬ 
son scattering and rie is the number density of electron. 


2.1 Ionizing photons from mini-quasars 

In this section, we derive the number of ionizing photons 
from the central accreting IMBHs, shining as mini-quasars. 
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Cosmic reionization of hydrogen and helium 3 


The ionizing radiation from accreting BHs primarily de¬ 
pends on the following three quantities: the redshift evo¬ 
lution of the mass density of BHs, the emission spectrum of 
the accreting BH and the duty cycle of BHs. 

Our calculation for the mass density of central BHs as 
a function of reds hift closely follows the simple model of 
IWang et al.l (l2010l l which is based on the PS mass function 
of dark matter halos. Here we only briefly describe the key 
features of this model, and refer the reader to that paper 
for more details. We assume that, at the center of each dark 
matter halo, there exists a BH with a mass (Mbh) that 
is correlated with the host halo mass. Then the mass den¬ 
sity of BHs could be estimated simply by calculating the 
number density of halos as a function of redshift. Here we 
ad opt the modified versio n of PS mass function proposed 
bv ISheth fc TormerJ (Il999l l to estimate the number of dark 
matter halos, which is 


nsT{M, z) AM 


^ 1 2 ui pm 

\ TT M 


1 -k 




a‘^{M,z) 


p-i 


Sc 


aiSc 

dlncr”'^ 


{M,z) 


AM 


{M,z) 

dM(3) 


where A — 0.3222, ai = 0.707, p = 0.3, Sc = 1.686 and pm 
is the current mean density of the Universe. The deviation 
of the linear density field (t(M, z) is given by 


CP{k)W^{k,M)Ak, (4) 

where W{k, M) is the top-hat filter, and D{z) is the growth 
factor of linear perturbations. 

The mass of the central BH (Mbh) is assumed to be 
correlated with the the stellar velocity dispersion cr* through 
the relation 


logMBH = + <)>log(CT*/CTo) , (5) 

where (Jq = 200kms~^. Basing on the improved 

self-adaptive mer g er tree model of dark matter halos 
llYuan et al.ll2008h . IWang et al.l ll2010ll explored the Mbh — 
cr* relation at high redshifts with the following BH forma¬ 
tion scenario: (i) each halo is randomly assigned a seed 
BH; (ii) seed BHs as remnants of Pop III stars have a 
slightly top-heavy initial mass function in the mass range 
200 Mq < Mbh < lOOOM© 0; (iii) the mass growth of BHs 
is considered to be dependent on gas accretion during halo 
mergers, and the mass accretion rate is Mbh = MBn/tef, 
where tef is the Salpeter e-folding time; (iv) the free param¬ 
eters, such as the radiative efficiency and Eddington ratio, 
are constrained to make sure the final BH mass agrees with 
the local Mbh — cr* relation. They found that the Mbh — cr* 
relation with parameters 9 = 6.67 and (j) = 2.79 can pro¬ 
vide a good fit to the numerical results at 2 : > 6. The stellar 
velocity dis persion cr* is related to the circular velocity of 
halos Vc by (lFerraresell2002l : IPlzzella et akll^OOSl l 

log Vc = 0.84 log a* -b 0.55. (6) 


^ Although recent studies imply a smaller characteristic mass of 
Pop III stars, it is yet inconclusive how this scena rio would affect 
the g rowth of MBHs. Besides, several authors fe.g. lVolonteri et al.l 
l2003h also argued that the evolution of MBHs is not very sensitive 
to the initial mass of seed BHs. Thus, we prefer a traditional seed 
BH model in this work. 


The relation betwee n the circular velocity 14 and the halo 
mass M is given by dBarkana fc Loebll200ll l 


K = 23.4 


with 


M 


10»/i-iM© 


1/3 


fin 


187r2 


1/6 


1 + z 
10 


1/2 


kms 

(7) 


Ac = 187r^ -b 82d - (8) 


and d = flin(-t) — 1. 

In Fig. [ 1 ] we show the BH mass densities as a function 
of redshift for IMBHs and massive BHs (MBHs; > 10® M©), 
respectively. We henceforth refer to accreting MBHs as 
“quasars”. For comparison, we also plot the strong upper 
limit on the IMBH densi ty at 2 > 6 of Pimbh < 3.8 x 
10'* M© Mpc“® derived bv ISalvaterra et ^ |2TO^, which is 
based on the measurement of the unresolved fraction of the 
SXRB. As can be seen from Fig. [T] the result is well below 
the constraints from the SXRB. 

Although the exact emission spectrum of mini-quasars 
remains poorly understood. Its typical spectrum Is expected 
to be harder than that of quasars. As suggested by ob¬ 
servations of ult raluminous X-ray sources (possibly associ¬ 
ated to IMBHs: iMiller et al.l 120031 In nearby galaxies, the 
mini-quasar spectrum could be approxlmatively described 
by two continuum comp onents, includ i ng b oth a “mul¬ 
ticolor disk blackbody” llMltsuda et al.l Il984h component 
and a non-thermal power-law component Le oc E~°‘ with 
Q sa 1. To simplify the calculation, we follow previous studie s 
dKuhlen fc Madaull2005l : lDiikstrall2006l : IZaroubi et al.l[200^ 

and assume that the mini-quasar spectrum has a simple 
power-law form. Specifically, two types of power-law spectral 
energy distributions (SEDs) are considered: the High Energy 
(HE) case with the energy range 200 eV < E < 100 keV 
and the Low Energy (LE) case with a low-energy cutoff of 
10.4 eV. Note that, the HE case assumes that all the ionizing 
UV photons are absorbed in the close vicinity of the source. 
As a comparison, we also note that the typical spectral in¬ 
dex of quasars is meas ured to be a = 1.7 (e.g. iTelfer et al.l 
I2OO2I: iTozzi et al.lf2006ll . 

Then the number of hydrogen (helium) ionizing photons 
from accreting IMBHs can be evaluated by 


dt 


L 


F{E)^ , 


JF{E)AE 


eLP, 


(9) 


where e Is the Eddington ratio, P represents the duty cycle 
and L Is the total luminosity emitted by accreting IMBHs 
per unit of comoving volume: 

L = j A'(M(Mbh, 2 ), 2 )LBdd(MBH)dMBH 
J 200 

= 1.26 X io38 Pim^h( 2:) ergs~*Mpc~®, (10) 

where I/Bdd is the Eddington luminosity, LEdd(MBH) = 
1.26 X 10 ®®(Mbh/M©) ergs“*, and P is the duty cycle 
of BHs, i. e. the probability that a BH is in the active 
state. A number of studies suggest that the d uty cycle in¬ 
creas e s rapidly with i n creasing reds hift (e.g. IWang et al.l 
I2OO6I : IShen et al.l l2007l : iTanaka 2014h. From t he measure- 
ments of clustering of quasars, Shankar et al.l (l2010l l sug¬ 
gested that P ~ 0.2, 0.5 and 0.9 at 2 = 3.1, 4.5 and 6, respec- 
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tively. To match these results, we therefore adopt a parame¬ 
terized form of P[z) = 0.03 + 0.97(z/2:tran)^/[l -f (z/ztran)^] 
with P = 2 and ztran = 4. To be optimistic, we also assume 
that all the ionizing photons produced by accreting BHs es¬ 
cape the host halo, i. e. /esc.BH = 1- Based on the above 
assumptions and descriptions, consequently, it is straight¬ 
forward to calculate the number of ionizing photons emit¬ 
ted at a particular energy range [iSmin, B'max] from accreting 
IMBHs per unit time. 


2.2 Ionizing photons from stars 

In order to determine the number of ionizing pho¬ 
tons produced by stars, the high-redshift star forma- 
tion rate (SFR) density is estima t ed ac cording to Q 
(|2008[1 . in which iHonkins fc BeacomI ll200(il Vs compilation 
on the SFR den s ity wa s updated with the ob s ervations by 
iGiavalis co et al. ll2004 l. Bunker e t al.l Ji ooll. lOuchi et al.l 


( 2004ll . iReddv et al.l 1 20081 1 and iBouwens et al l ( 2008l l. Bas- 
ing on these observational d ata, O (l2008ll obt a ined a 
piecewise-linear fitting, as in iHookins &: BeacomI ll2006ll . 
which is 


logp«(a) = a-I-folog(l-I- 2 ), (11) 


where 

['(-1.70,3.30), 2 <0.993 

(a, 6 ) = i (-0.727, 0.0549), 0.993 < 2 < 3.80 , (12) 

[(2.35,-4.46), 2 > 3.80 

and p* is the SFR in units of M© yr”*^ Mpc“®. Then the 
number of ionizing photons is calculated by 

— n-,./escP*(2), (13) 

where /esc represents the escape fraction of ionizing pho¬ 
tons from stars and n.y gives the number of photons emitted 
per unit mass of stars. Unfortunately, the escape fraction of 
stellar ionizing photons from high-redshift galaxies remains 
poorly constrained, making this parameter one of the pri¬ 
mary goals of many theo retical and observational studies. 
A number of studies (e. g. Ilwata et al.ll200^ : IVanzella et al.l 
l2nifll : iNestor et al.ll20llh have measured the escaping ioniz¬ 
ing photons at 2 ~ 3 — 4, suggesting that the escape fraction 
is between 0.05 and 0.2. Although the samples are small, 
almost all the average escape fractions of ionizing photons 
from the observed massive galaxies are less than 0.2. By 
integrating over appropriate energy ranges, the parameter 
n-^ gives (8.05, 2.62, 0.01) x 1 0^° M~^ for HII, Hell , Bell i, 
respectively, as suggested bv IChoudhurv fc Ferraral ll2005ll . 

For comparison, we show how the ratio 
evolves as a function of redshift for three different /esc in 
Fig. [21 Here we consider only the LE case for the IMBH 
spectrum, as for the HE case, the hydrogen ionizing photons 
from mini-quasars are simply zero. It is also worth stressing 
that all the ionizing photons from mini-quasars are assumed 
to be available for ionizing the IGM. As shown in Fig. (21 the 
ratio peaks at 2 : ~ 13 for all three models, decreases rapidly 
at lower redshifts and becomes negligible at 2 < 5. Only 
for the model with /esc < 0 . 2 , mini-quasars have a higher 
ionizing flux contribution than that of stars at 2 > 6 . 


3 RESULTS 


We first consider the contribution of mini-quasars, i. e. ac¬ 
creting IMBHs alone to the reionization of cosmic hydrogen 
and helium. Fig. [3] shows the evolution histories of the fill¬ 
ing factor of the ionized regions Qhii (solid line) and (Jneiii 
(dashed line) due to the accreting IMBHs alone. Note that 
only the LE case for the IMBH emission is shown here, since 
accreting IMBHs only provide a negligible contribution to 
hydrogen reionization for the HE case as mentioned above. 
As can be seen, the accreting IMBHs can only contribute 
~ 20% of hydrogen ionizing photons at 2 ~ 6 . The differ¬ 
ence between the evolution of Queiii and that of Qhii is 
very subtle, because of the assumption that the He HI ion¬ 
ization front never overtakes the hydrogen ionization front, 
i. e. Qhii ^ QHeiii. This result implies that, in an opti¬ 
mistic case (LE case), the contribution of mini-quasars to 
hydrogen reionization is not negligible at 2 ~ 6 , but still 
insufficient to provide the required num ber of ionizing pho¬ 
tons, consistent with previous studies dMadau e t al.l 2004) : 
iDiikstra et al.l[2004l : IWang et al.ll201ol : iMcQuinn 2012lL 


As a reference, the contribution from accreting MBHs 
(“quasars”) is also considered. Note that we set the emis¬ 
sion spectrum of quasars with a spectral index of a = 1.7 
and ionizing UV and X-ray photons (i. e. energy range of 
10.4eV < E < 100 keV). As shown in Fig. O only at rel¬ 
atively low redshifts (2 < 6 ) does the contribution from 
MBHs play a relevant role in reionizing the IGM. When 
considering the contributions from both IMBHs and MBHs, 
we find that helium reionization ends at 2 ~ 3, consis- 


Zhene et al. 

l2004l:lMcOuinrj2009l:lFurlanetto & DixorJ201C 

Shull et al. 

2 OIOI: Becker et al.l 201 ll: Worseck et al.ll 2011 


middle of helium reionization occurs at 2 ~ 5 and hydro¬ 
gen reionization ends at 2 « 4. It is worth stressing that, for 
simplicity, we have assumed the same Mbh — cr« relation at 
all redshifts. This almost certainly constitutes an underesti¬ 
mate for the contribution of accreting BHs at low redshifts 
and delays the completion of helium reionization. 


We then consider the combined contribution from ac¬ 
creting BHs and stars to the reionization of hydrogen and 
helium. Here, for the ionizing emission from accreting BHs, 
both IMBHs and MBHs are included. Given the uncer¬ 
tainty in the escape fraction of ionizing photons from stars, 
we present the results of modeling the reionization histo¬ 
ries of hydrogen and helium for three different choices of 
the escape fraction of ionizing photons from stars (/esc = 
0.05; 0.2; 0.5). Fig. |4]shows the evolution of both Qhii and 
Queiii- The upper panel refers to the HE case for the IMBH 
emission, while the lower panel is for the LE case. It is clear 
that the Universe is reionized earlier for the LE case. As seen 
in both panels, the escape fraction /esc of ionizing photons 
from stars does not have strong influence on the evolution 
of Queiii, because stars provide a negligible contribution to 
the helium ionizing photons. 


In the HE case, helium reionization ends at 2 « 2.7, 
consistent with observations, while in the LE case the end 
redshift of helium reionization 2 « 3.5 is much earlier. Mod¬ 
els with /esc = 0.2 and /esc = 0.5 satisfy the observational 
constraints of the Gunn-Peterson trough (i. e. hydrogen 
reionization at 6 < 2 reion Pz 7). However, when considering 
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the constraints on the CMB optical depth, only the model 
with /esc = 0.5 is consistent with the WMAP-7 result of 
Te = 0.088 ± 0.015, as shown in Fig. O 

In the LE case, only the model with /esc = 0.2 is 
consistent with observations of the Gunn-Peterson trough, 
whereas hydrogen reionization in the other two models is 
completed either too early {z > 7) or too late {z < 5.5). 
However, although the model with /esc = 0.2 is able to 
reionize the Universe at 2 6, it is not able to satisfy the 

constraints on the CMB optical depth simultaneously. This 
result imply that, to simultaneously satisfy the joint con¬ 
straints from the Gunn-Peterson test and the CMB optical 
depth, an escape fraction /esc within the range from 0.2 to 
0.5 is needed in this case. 

So far, we have analyzed the evolution of the reion¬ 
ization histories for different escape fractions, assuming a 
power-law spectral index of q = 1 for the ionizing emis¬ 
sion of accreting IMBHs. For comparison, we also explore 
the possibility of a quasar-like spectral index of a = 1.7 for 
the emission spectrum of IMBHs. The results are shown in 
Fig. [6]and Fig.[7l As seen in Fig. [6l helium reionization ends 
much earlier at z > 3.5 for both HE and LE cases. For the 
HE case, the results turn out to be similar to those found 
in the a = 1 scenario, except that helium reionization ends 
much earlier at a > 3.5. On the other hand, for the LE case, 
all three models with different /esc can produce enough ion¬ 
izing photons to obtain the CMB optical depth, as shown 
in Fig. [7] In this extreme case, only an escape fraction of 
/esc < 0.05 is large enough. These results show that, if high- 
redshift star-forming galaxies have an escape fraction larger 
than 0.05, then accreting IMBHs with a quasar-like spec¬ 
trum would produce too many ionizing photons to be more 
or less in tension with observations. 

In the most optimistic model, mini-quasars are able to 
reionize more than cnj 20 % of the cosmic helium by redshift 
z ~ 6, and the ionized fraction of cosmic helium rapidly 
climbs to more than 50% by redshift z 5 when consider¬ 
ing both IMBHs and MBHs. These results may imply that, 
in addition to the constraint from the SXRB, better mea¬ 
surements of the evolution history of helium reionization, 
especially at high redshifts, would also be helpful in placing 
limits on the growth of IMBHs in the early Universe. If suc¬ 
cessful, this would not only represent a breakthrough in our 
understanding of the epoch of reionization, but it would also 
shed some light on the puzzles concerning the formation of 
SMBHs at 2 ~ 7. 


4 DISCUSSION AND CONCLUSIONS 

In this paper, basing on recent observational results, we have 
studied the reionization history of cosmic hydrogen and he¬ 
lium by both mini-quasars and stars. The contribution from 
quasars is also included. In order to take into account the 
contribution from (mini-)quasars, the Mbh — o'* relation at 
z > 6, which is from the simulation on the merger tree of 
dark matter halos, is applied. Several types of power-law 
spectra for the mini-quasar emission are considered. Three 
different escape fractions /esc of ionizing photons from stars 
are tested. By simultaneously considering the constraints 
from the CMB optical depth and observations of the Gunn- 
Peterson trough in the spectra of high-redshift quasars, we 


find the following results for models with a power-law IMBH 
spectrum (spectral index a = 1). 

(i) For mini-quasars with no UV ionizing photons 
(200 eV < E < 100 keV), the hydrogen ionizing photons 
emitted by mini-quasars are negligible. Helium reionization 
occurs at 2 « 2.7, consistent with observational results. Only 
the model with /esc ~ 0.5 is simultaneously consistent with 
the constraints from the Gunn-Peterson test and the CMB 
optical depth. 

(ii) When considering mini-quasars with UV ionizing 
photons (10.4eV < E < 100keV), both hydrogen and he¬ 
lium reionization are completed earlier. In this case, an es¬ 
cape fraction in the range 0.2 < /esc < 0.5 is able to sat¬ 
isfy the joint constraints. Note that, in this optimistic case, 
mini-quasars are able to contribute ~ 20% of h ydrogen ioniz¬ 
ing photons, consistent with previous studies dMadau et al.l 
12004 IWang et al.ll2010ll . 


We also consider the possibility that the mini-quasar has a 
quasar-like spectral index of q = 1.7. In this scenario, mini¬ 
quasars significantly over-produce the helium ionizing pho¬ 
tons, leading to an early completion of helium reionization 
at z > 4:, which is inconsistent with observations. 

There are some caveats which could significantly affect 
our expectations about the contribution to cosmic reioniza¬ 
tion from mini-quasars. First, we have assumed that all the 
ionizing photons from mini-quasars are available for ionizing 
the IGM. This is because the high production rate of ion¬ 
izing photons from mini-quasars woul d allow the surround - 
ing gas t o be ionized ve r y qui ckly fe.g. IWhalen et ^|2004| L 
However. iBenson et al.l ( 2013l l argued that the escape frac¬ 
tion of ionizing photons from the early quasars could be 
only /esc.BH ~ 0.3. If this is the case, the contribution from 
mini-quasars would be significantly reduced. Another im¬ 
portant consideration is the assumption on the duty cycle 
of mini-quasars at high redshifts, i. e. P{z > 6) ~ 0.9 — 1. If 
the duty cycle of mini-quasars at high redshifts is similar to 
that of present quasars {P{z ~ 0) ~ 0.1), the contribution 
of mini-quasars to cosmic reionization would become negli¬ 
gible. Besides, on the assumption of the spectrum of mini¬ 
quasars, instead of a simple powe r-law, one can also adopt 
a multicomponent spectrum fe.g. iTanaka et alll2012h . Also 
note that we have ignored the secondary ionizations from 
the mini-quasar emission in our model for simplicity. We 
caution that, until most of these properties of mini-quasars 
have been better understood, a stringent constraint on the 
contribution of mini-quasars can not be obtained. In turn, 
it is also reasonable to expect that, improved measurements 
of helium reionization, especially at high redshifts, coupled 
with the constraint from the SXRB, could be used to bet¬ 
ter constrain these properties of mini-quasars, which would 
make a relevant contribution to the solution of the puzzles 
related to the formation of SMBHs at high redshifts. 

Moreover, it is commonly believed that the feedback 
from the accretion onto central BHs could play an impor¬ 
tant role. For example, the radiation from accreting BHs 
can limit the gas supply by radiation pressures and photon- 
heati ng, leading to the self-regulation of t he BH growth 
(e.g. Jgon__etaL 20121 : [^ark fc Ricottil [201^ 1. Additionally, 
iTanaka et al.l ( 2012l ~) suggest that the IGM heating by mini¬ 
quasars can also suppress the formation and growth of cen¬ 
tral BHs in low-mass halos. On the other hand, the star for- 
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mation in low-mass galaxies may be affect ed by both local 
and global heating fro m mini-quasars (e.g. iRipamonti et al] 
l2008l : [jeon et al1l2012li . This would thus influence the rela¬ 
tive contributions of mini-quasars and stars to cosmic reion¬ 
ization. However, this additional complication is beyond the 
scope of this paper, and will be addressed in our future work. 
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Figure 1. Evolution of the mass density of central BHs as a function of redshift. The solid line and short-dashed line represent the 
mass density of I MBH s and MBHs, respectively. The dashed line is the constraints on the BH mass density from the SXRB derived by 
ISalvaterra et all ll2005h . 



Figure 2. Redshift evolution of the ratio between the IMBH and star contributions to the hydrogen ionizing photons, assuming the LE 
SED (10.4 eV < E < 100 keV) for the IMBH emission. The escape fractions of ionizing photons from stars are taken to be /esc = 0.05, 
0.2 and 0.5 for the solid, dashed and short-dashed lines, respectively. 
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Figure 3. Redshift evolution of the volume filling factor of Qhii QHeiii due to the ionizing photons from IMBHs and MBHs, 
assuming the LE SED (10.4 eV < E < 100 keV) for the IMBH emission. 
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Figure 4. Redshift evolution of the volume filling factor of Qhii und QHeiii for three different escape fractions /esc of ionizing photons 
from stars. Here, both accreting BHs and stars are included. The IMBH spectrum has a spectral index of a = 1. The upper panel refers 
to the HE case for the BH emission, while the lower panel refers to the LE case. 
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Figure 5. Redshift evolution of the Thomson scattering optical depth for three different escape fractions /esc. The upper and lower panels 
refer to the HE and LE cases for the IMBH emission, respectively. The black bar represents the WMAP-7 result, i.e., Te = 0.088 di 0.015. 



Figure 6. Redshift evolution of the volume filling factor of Qhii and QHeiii> assuming the mini-quasar emission with a quasar-like 
spectrum (q = 1.7). The upper and lower panels refer to the HE and LE cases for the IMBH emission, respectively. Here, Hell 
reionization completes much earlier at z > 3.5 for both HE and LE cases. 
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Figure 7. Redshift evolution of the Thomson scattering optical depth, assuming the mini-quasar emission with a quasar-like spectrum 
(a = 1.7). The upper and lower panels refer to the HE and LE cases for the IMBH emission, respectively. The black bar represents the 
WMAP-7 result, i.e., = 0.088 ib 0.015. 
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